Anomalous Suppression of Valley Splittings in Lead Salt Nanocrystals without 

Inversion Center 



A.N. Poddubny,! M.O. Nestoklon/ S.V. Goupalov^'^ 
^ loffe Physical- Technical Institute, 26 Polytekhnicheskaya St., 194021 St. Petersburg, Russia 
^ Department of Physics, Jackson State University, Jackson, Mississippi 39217, USA 

Atomistic sp^d^s* tight-binding theory of PbSe and PbS nanocrystals is developed. It is demon- 
strated, that the valley splittings of confined electrons and holes strongly and peculiarly depend on 
the geometry of a nanocrystal. When the nanocrystal lacks a microscopic center of inversion and has 
Td symmetry, the splitting is strongly suppressed as compared to the more symmetric nanocrystals 
with Oh symmetry, having an inversion center. 
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Interest in almost spherical nanocrystals (NCs) made 
of lead chalcogenides (PbSe, PbS) has recently exploded 
due to their utility in fundamental studies of quantum 
confinement effects and enabling potential for technolog- 
ical applications in photovoltaics [IHl]- These applica- 
tions impose requirements of reproducibility of the de- 
vices, ability to control their properties and necessity 
to understand mechanisms behind physical processes un- 
derlying their work. It has been noted [5| [6J that lead 
salt NCs are very peculiar compared to quantum dots 
(QDs) of III-V and II- VI compound semiconductors be- 
cause lead chalcogenides have band extrema in the four 
inequivalent L-points of the Brillouin zone and such ef- 
fects as confinement-induced valley-mixing and effective 
mass anisotropy should be considered to fully account 
for the properties of lead salt NCs. It also was suggested 
that optical properties of lead salt NC QDs can be sen- 
sitive to a particular arrangement of atoms within the 
QD, as the latter determines the overall symmetry of the 
structure [71 [8]. 

In this work we study how all these effects influence the 
valley-orbit and spin-orbit splittings of one-particle en- 
ergy levels of electrons and holes confined in lead salt NCs 
using the atomistic tight-binding (TB) approach. We 
found that these splittings are very sensitive to the par- 
ticular arrangement of atoms within the almost spherical 
NC. In particular, we considered NCs of almost spheri- 
cal shape centered on an anion or cation atom, serving 
as a center of inversion, along with NCs having no inver- 
sion symmetry. We found that in NCs without a center of 
inversion the valley-orbit and spin-orbit splittings of elec- 
tron energy levels are strongly suppressed. This effect is 
quite unusual because typically a higher symmetry of a 
physical system implies a higher degeneracy of its energy 
levels, while in our case the suppression of the splittings 
occurs in NCs having lower symmetry. Nevertheless, we 
were able to explain this puzzling behavior using mathe- 
matical apparatus of the group theory. 

Lead chalcogenides (PbSe, PbS) are semiconductor 
compounds with a rocksalt crystal lattice and a narrow 
and direct band gap [9 . The extrema of both the conduc- 



tion and valence bands are located at the four L-points 
of the Brillouin zone: 

fcl,2 = f (1, ±1, ±1) , fc3,4 = f (-1, ±1, Tl) , (1) 

where a is the lattice constant. The empirical TB method 
is an efficient tool to model electronic properties of large- 
scale nanostructures ^OJ . Success of the TB parametriza- 
tion depends on the choice of basis functions and on 
the accuracy of the fit of the bulk band structure. The 
simplest TB parametrizations of lead chalcogenides are 
based on the basis set of the three p orbitals playing ma- 
jor role in the formation of the valence and conduction 
band states [TT1[T2]. More quantitatively accurate mod- 
els include also 5* and d orbitals ^ ITSHTS] . However, 
no attempts (with the only exception of Ref. [12]) have 
been made to fit the actual effective masses of the elec- 
trons and holes near the L-points. On the other hand, 
the second-nearest neighbors model of Ref. [12] fails 
to reproduce the bulk dispersion for wavevectors far from 
the L points [6 . Consequently, even the most advanced 
existing TB parametrizations of lead chalcogenides [5 are 
not suitable [16 for an adequate description of the NCs. 

We have performed an independent atomistic sp^d^s* 
TB parametrization of the electron energy dispersion in 
bulk PbSe and PbS by fitting the spectra calculated by 
the state-of-the art GW technique of Ref. [17\ The 
goal values for the carrier effective masses near the L- 
points were set to the experimental values [18 : rn^^l^ = 
0.070 mo, mlj 0.040 mo, m^J 0.068 mo, m^J 
0.034 mo for PbSe and rrf^f^ = 0.105 mo, Tn^Y = 
0.080 mo, mlj = 0.105 mo, m^J = 0.075 mo for PbS 
(mo is the free electron mass), as even the modern ab 
initio approach [17 does not satisfactory reproduce the 
effective masses. 

The TB parameters we obtained are listed in Table [Tj 
The resulting effective masses ruc^i = 0.068 mo, mc,t = 
0.041 mo, rrih^i = 0.069 mo, rrih^t = 0.039 mo for 
PbSe and rric^i = 0.098 mo, mc,t = 0.079 mo, m^^^ = 
0.104 mo, rrih^t = 0.074 mo for PbS are quite close to the 
experimental values. The spin-orbit coupling constants 
of p orbitals at Pb, Se, and S were not changed during the 
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FIG. 1: Central parts of the three types of nanocrystals. 



fitting procedure and were taken from Refs. [21] and [22] 
for Pb and for the anions, respectively. 

In our study we will consider the three types of NC 
geometries illustrated in Fig. [ij For the structures 
shown in Figs.[l^(b) the center of the spherical NC is on a 
cation (anion) atom while for the structure in Fig. ^ the 
center of the sphere lies halfway between a cation and 
an anion on a line parallel to the [111] direction. The 
non-stoichiometric QDs of the types (a) and (b) both 
have centers of inversion and are characterized by the 
cube symmetry group Oh- The stoichiometric QD of the 
type (c) has no inversion center and is characterized by 
the tetrahedron symmetry group T^. Note that the QDs 
cannot be perfectly spherical due to the discretness and 
lower point symmetry of the underlying crystal lattice. 
In our work the QDs are formed by all the atoms within 
a certain distance from the center of the NC. It is conve- 
nient to measure this distance with a dimensionless inte- 
ger number. Thus, we define the "number of shells" as 



TABLE I: TB parameters for PbSe and PbS. The transfer 
integrals are meausured in eV and given in the Slater-Koster 
notations . The spin-orbit splittings are defined according 
to Ref. EOl. 
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FIG. 2: (Color online) Energy levels in PbSe (a,b,c) and PbS 
(a^,b^,c^) NCs with diameter D ^ 4.9 nm and D ^ 4.6 nm, 
respectively. Panels (a-c) and (a^-c^) correspond to cation- 
centered NCs, anion-centered NCs and NCs without an inver- 
sion center, respectively (see Fig. [Tj). States of the odd and 
even parity are marked by the long thick blue and short thick 
red lines, respectively. The thin black lines correspond to the 
states without a certain parity. The degeneracies of the low- 
est electron and hole confined states are indicated near each 
line. 



the number of atomic layers within the distance from the 
center of the QD to its surface along the [100] direction. 

Contrary to covalent semiconductors like Si [5 , lead 
chalcogenides are characterized by the strongly ionic 
atomic bonds making them relatively insensitive to the 
surface chemistry [16j. In particular, no surface states 
appear in the fundamental band gap of non-passivated 
lead chalcogenide NCs. Therefore, we have not passi- 
vated the surface atoms in our TB modeling. In real 
QDs of the types (a) and (b) (cf. Fig.[T]) such passivation 
is necessary to compensate for the surface charge [23] . 
The actual structure of the NCs should depend on the 
details of the synthesis procedure and can be determined 
with the help of the nuclear magnetic resonance [24 and 
X-ray diffraction [25 techniques. 

The calculated energy levels of confined carriers for 
PbSe and PbS NCs of the diameter D ^ 5 nm (cor- 
responding to 9 shells) are shown in Fig. |2] For each 
material three panels (a-c and a'-c') correspond to the 
three possible NC geometries illustrated in Fig. [TJ The 
band gap in both cases agrees well with the results of 
Ref. [5.. 

All the states can be divided into distinct groups char- 
acterized by a certain parity. For NCs with a center of 
inversion, each state automatically has a certain parity. 
Indeed, in bulk lead chalcogenides the lowest electron 
state in the conduction band has the Lq symmetry [17], 
i.e. it is odd with respect to the inversion symmetry 
operation when the center of inversion is chosen on the 
cation atom. The uppermost electron state in the va- 
lence band has the opposite parity. For QDs without an 
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inversion center one can define approximate projectors to 
the even and odd states. We have attributed a certain 
parity to the states, for which the squared mean values 
of the projectors differed more than in three times. 

Energy splittings within each multiplet characterized 
by a certain parity are clearly seen in Fig.|2]and can be ex- 
plained by the confinement-induced inter- valley coupling 
and the carrier effective mass anisotropy. The impor- 
tance of these two effects for lead chalcogenide NCs has 
been emphasized in Ref. [6 . However, the dependence 
of the splittings on the NC geometry clearly manifested 
m Fig. [2] has never been reported. 

A striking feature of Fig. |2] is the suppression of the 
energy splittings for the type (c) NCs lacking a center 
of inversion. The splittings are quite small and cannot 
be distinguished within the energy scale of Fig [2] On 
the contrary, for QDs with an inversion center [panels 
(a), (b), (a^), (b^) of Fig. [2], the ground-state multiplets 
for both electrons and holes have well-pronounced struc- 
tures with substantial splittings even for QD diameters 
as large as 4.9 nm. This observation refers to both the 
conduction and valence band electron states. The effect 
is more pronounced for the PbS QDs than for the PbSe 
QDs, which can be related to the more isotropic effective 
masses of the band extrema in bulk PbS. 

To elucidate this puzzling behavior, we have analyzed 
the dependence of the splittings on the NC diameter. 
For simplicity, we restrict our consideration by the elec- 
tron and hole ground states. Within the effective mass 
approximation, the ground state of confined carriers is 
fourfold degenerate with respect to the valley index and 
twofold degenerate with respect to the spin projection, 
i.e. the total degeneracy is eightfold. If we neglect 
the spin and consider valley-orbit interaction only, the 
ground state is split into a state of Ai symmetry (sin- 
glet), and a state of F2 symmetry (triplet), as sketched 
in insets of Figs. 3j4 When the spin degree of free- 
dom is taken into account then both the singlet and the 
triplet states acquire extra degeneracy. This degeneracy 
is partly lifted, as the six-fold degenerate state corre- 
sponding to the triplet is split by the spin-orbit inter- 
action into a two-fold degenerate state of E2 symmetry 
and a four- fold degenerate state of G' symmetry [26] . As 
a result, the carrier ground- state level is split into the 
three multiplets: the two doublets (of E[ and E2 sym- 
metry, respectively) and the four-fold degenerate state 
of symmetry. As far as the symmetry with respect 
to inversion is not concerned, the symmetry groups 
and Oh are equivalent. Therefore, this symmetry anal- 
ysis applies to all types of NC geometries presented in 
Fig.[l] 

Figures |3] and [4] show the energies of the resulting con- 
duction (valence) band multiplets in PbS NCs as func- 
tions of the NC diameter. The panels (a)-(c) correspond 
to the three NC geometries considered throughout the 
paper (see Fig.[l]). The energies of the states are counted 
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FIG. 3: (Color online) Energies of levels belonging to the 
ground state multiplet of the conduction band electron in PbS 
NCs as functions of NC diameter. Panels (a), (b), (c) corre- 
spond to Pb-centered NCs, S-centered NCs and NCs without 
inversion center, see Fig. [l] Squares, triangles and circles 
correspond to the states with the symmetry G' , E2 and E[, 
respectively, see the level splitting scheme in the inset. 



from the averaged value {Ee[ + Ee'^ + 2Eg')/4. The 
splittings strongly oscillate with the number of shells 
A/" in a NC. Such oscillations are typical for the val- 
ley splittings in various semiconductor structures. Sim- 
ilar behavior has been reported for SiGe/Si[27l EH] and 
GaSb/AlAs[29l[30] quantum wells and Si NCs [ST. 

Comparison of panels (a) and (b) of Figs. |3|4| on one 
hand with the panels (c) of Figs. 3|4 



on the other hand 

clearly shows that the suppression of valley splittings in 
NCs without a center of inversion is a general feature per- 
sistent in a wide range of NC sizes. Comparison of Fig. [3] 
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FIG. 4: (Color online) Same as Fig. [s] but for the valence- 
band ground state. 
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with Fig.[4]enables one to conclude that the spin-orbit in- 
teraction is much stronger for conduction-band electrons 
than for valence-band holes. Indeed, the energies of the 
hole states with the symmetry G' and E2 in Fig. |4] are 
almost the same, while, for conduction-band electrons, 
all the splittings in Fig. [3] are of the same order. 

Let us explain the anomalous suppression of the valley 
splittings in lead salt NCs without an inversion center. 
We want to account for the inter-valley coupling in the 
lowest non- vanishing approximation. To this end we con- 
sider electronic states originating from the four inequiva- 
lent L valleys of a bulk semiconductor and neglect the A;-p 
mixing of conduction and valence band states [9 . Then 
the wave function of the confined electron state in the j- 
th L- valley can be written as {r\Lj) = e^^^'^Uj{r)^j{r) , 
where ^j{r) is the smooth envelope function, Uj{r) is 
the periodic Bloch amplitude for the bulk state in the 
j-th valley, and spin indices are omitted. We will fur- 
ther assume that the bulk material has isotropic effective 
masses of the band extrema in L points. In this approxi- 
mation the envelopes ^j{r) are invariant under rotations. 
The confinement-induced inter- valley coupling can be de- 
scribed by the following matrix element: 

Ij^k = {Lj\HQj,\Lk), j,/c = 1...4, j^k, (2) 

where i^QD is the microscopic QD Hamiltonian. Then it 
follows that the integral Ij^k vanishes when the QD lacks 
inversion symmetry, i.e. belongs to the type (c). To show 
this let us rewrite L k as 



3,k 
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-ik 



^u;{r)H,uike''''^Uk{r)dr 



where the integral in the right-hand side is over a unit 
cell and contains the Hamiltonian of a bulk material while 
the summation runs over all the cation (or anion) sites 
within the QD. It is this summation that is sensitive to 
the arrangement of atoms within the QD. For the type (c) 
geometry the sum is exactly zero. This cancellation takes 
place independently of the radius of the QD and is fully 
determined by the symmetry. To see this one can use the 
following well known fact [26] . If a given function describ- 
ing some crystalline physical system transforms accord- 
ing to a certain representation of the system's symmetry 
group, then the sum of this function over the lattice sites 
belonging to the system may be different from zero if 
and only if the decomposition of this representation into 
irreducible ones contains the identity representation. 

In our case one can distinguish three linearly indepen- 
dent functions exp[i{kk — kj)Rn] which may be chosen 
as shown in the first column of Table |TI[ Table [TT] gives 
the values of these exponent functions when Rn sweeps 
the coordinates of the anion atoms shown in Fig. [ijc). 



These atoms may be obtained from one another by the 
rotations of the type (c) QD. The last column of Ta- 
ble |ll| indicates that the exponent functions transform 
according to the vector irreducible representation F2 of 
the group Td. This representation is different from the 
identity representation Ai. Therefore, for type (c) QDs 
Eq. (§ is zero. Table III skives the values of the same 



exponent functions when Rn sweeps the coordinates of 
the anion atoms shown in Fig.jlja). These atoms may be 
obtained from one another by the rotations of the type 
(a) QD. The last row of Table |III| shows that the sum 
of the exponent functions remains invariant under such 
rotations. More precisely, the exponent functions trans- 
form according to the direct sum of the two irreducible 
representations Af ®E~^ of the group Oh- Thus, for type 
(a) QDs Eq. ([3| is different from zero. 

TABLE II: Phase factors in Eq. ([3| for the coordinates of the 
anion atoms in Fig. [TJc). 
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TABLE III: Same as Table [TT| but for anion atoms in Fig.[T]^8 
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This consideration is no longer valid if the function 
^j{r) is anisotropic. This is the case of real lead salts 
NCs, as in bulk lead chalcogenides the longitudinal mass 
in the L valley is larger than the transverse one. Con- 
sequently, in real NCs lacking inversion center the valley 
splitting is not exactly zero but determined by the de- 
gree of the effective mass anisotropy in L valleys. This 
explains the fact that in PbS NCs the splitting is smaller 
than in PbSe ones, cf. panels (c) and (c^) of Fig. [2] 

In conclusion, we obtained a new set ofsp^d^s* TB pa- 
rameters for the bulk PbSe and PbS semiconductor com- 
pounds and calculated the electron and hole energy levels 
in NCs made of these materials. We demonstrated that 
the valley-orbit and spin-orbit splittings of the ground 
electron and hole energy levels are very sensitive to a 
particular arrangement of atoms in the NC and can be 
strongly suppressed for a certain geometry, when the NC 
lacks a center of inversion. 
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